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ABSTRACT 

In  [l]  and  [2]  identities  for  the  wave  equation  similar  to 
the  energy  identity  were  derived  by  showing  that  if  u  is  suf- 
ficiently smooth,  there  exist  linear  first  order  operators  Nu 
such  that  Nu  Du,  with  Q  =  S  /St   -  A,  can  be  written  as  the  di- 
vergence of  a  vector.   Thus  if  Q  u  =  0  a  certain  surface  integral 

in  (x, t)  space  vanishes  identically  and  this  yields  the  identity. 

2   2 
One  of  these  operators  is  2rtu  +  ( r  +t  )u,  +  2tu,  another  is 

ru  +  tu,  +  u  with  r  =    |x| .   For  the  familiar  energy  identity 
Nu  =  u,  . 

In  Part  I  of  this  report  these  identities  will  be  re- 
derived  for  three  space  variables  by  noting  that  certain  trans- 
formations, in  particular  the  Kelvin  transformation,  leave  the 
wave  operator  invariant  and  hence  the  classical  energy  identity 
can  be  transformed  into  other  identities. 

In  Part  II,  the  Kelvin  transformation  and  the  resulting 
identity  are  applied  to  incoming  and  outgoing  waves  as  defined 
by  Lax  and  Phillips  [5>^]- 

In  Part  III  the  main  theorem  of  the  first  part  is  used  to 
prove  the  following  result  in  geometrical  optics:   Suppose 
that  v;e  are  given  a  smooth,  star-shaped  perfectly  reflecting, 
three-dimensional  body  that  extends  to  infinity  and  that  a  high 
frequency  harmonic  source  of  light  illuminates  the  region  out-. 
side  the  body  in  such  a  way  that  no  shadow  is  cast.   The  field 
is  given  by  a  solution  of  a  boundary  value  problem  for  the  re- 
duced v.'ave  equation.   There  is  also  an  approximate  solution 
given  by  geometrical  optics.   The  theorem  states  that  these  tv/o 
are  asymptotically  equal  in  the  limit  of  infinite  frequency  for 
the  harmonic  source. 
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In  [1]  and  [2]  identities  for  the  wave  equation  similar 
to  the  energy  identity  were  derived  by  showing  that  if  u 

is  sufficiently  smooth,  there  exist  linear  first  order 

2    2 
operators  Nu  such  that  NuDu,  with  D  =  o  /ot  -  A,  can 

be  v;ritten  as  the  divergence  of  a  vector.   Thus  if  D  u  =  0 

a  certain  surface  integral  in  (x,t)  space  vanishes  identically 

and  this  yields  the  identity.   One  of  these  operators  is 

2   2 
2rtu  +  (r  +t  )u,  +  2tu,  another  is  ru  +  tu,  +  u  with 
r  t  r     t 

r  =  l3c|  .   For  the  familiar  energy  identity  Nu  =  u ,  . 

In  Part  I  of  this  report  these  identities  will  be 
rederived  for  three  space  variables  by  noting  that  certain 
trans format  ions J  in  particular  the  Kelvin  transformation, 
leave  the  wave  operator  invariant  and  hence  the  classical 
energy  identity  can  be  transformed  into  other  identities. 

In  Part  II,  the  Kelvin  transforma-uLon  and  the  resulting 
identity  are  applied  to  incoming  and  outgoing  waves  as 
defined  by  Lax  and  Phillips  [p.^J. 

In  Part  III  the  main  theorem  of  the  first  part  is  used 
to  prove  the  following  result  in  geometrical  optics:   Suppose 
that  we  are  given  a  smooth,  star-shaped  perfectly  reflecting, 
three-dimensional  body  that  extends  to  infinity  and  that  a 
high  frequency  harmonic  source  of  light  illuminates  the 
region  outside  the  body  in  such  a  way  that  no  shadov;  is 


-1.1  arisll 


fix 
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cast.   The  field  is  given  by  a  solution  of  a  boundary  value 
problem  for  the  reduced  wave  equation.   There  is  also  an 
approximate  solution  given  by  geometrical  optics.   The 
theorem  states  that  these  two  are  asymptotically  equal  in 
the  limit  of  infinite  frequency  for  the  harmonic  source. 


I.   Energy  Identities 

1.   The  classical  energy  identity. 

2 
Lemma  1.   If  u  has  second  derivatives  in  L  in  D  then 

2j(Ju^  Quldxi  =  J  {{n^^+{]^u\^)t^^  -   2u^  u^  x^)do- 
D  D 

where  t  , x  are  the  components  of  the  unit  normal  in  time 
and  space  and  u  is  the  derivative  in  the  direction  of  the 
outward  space  normal  to  D. 

For  future  reference  we  note: 

a)  On  that  part  of  D  for  which  t  is  constant  the  inte- 
grand  is  u,  +  |Vu|  . 

b)  On  that  part  of  D  which  is  independent  of  time 
the  integrand  is  -2^^  u  . 

c)  On  that  pai't  of  D  which  is  characteristic,  x  =  +  t 

•^  n   —  n 

the  integrand  is  (u^  +  (jVul   -  u„))  where  u_  is  the 

derivative  along  the  bicharacteristic  in  Dj  u  is  the  outward 


3io   orIT 


•       ^^ 


-'fi;t:  ff*.  ! 
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space  normal  derivative.   Thus  the  integrand  involves  only 

derivatives  in  the  surface. 

It  is  convenient  to  have  the  three-dimensional  wave 

operator  in  polar  coordinates  (r,9/+)  and  to  introduce 

w  =  ru  as  dependent  variable.   Setting  .A  =  sin"  Ob  /b'i>      + 
-1. 


sin"  Q^lsin  05/D9)/Sa,  we  find 


-1        1  -2 

O  u  =  r   Lw  =  r   (w^^  -  w^^^  -  r  Aw) 


By  substitution  in  the  energy  identity  one  can  find 

an  identity  for  w.  A  more  convenient  expression  involving 

_2 
only  derivatives  of  w  is  obtained  by  adding  to  u ,  0  u  =  r  w.Lv/ 

a  divergence  expression  which  vanishes  identically.   In 

-1  2        -1  2 
this  case  we  add  (r  w  )  ,  -  (r  ^  ),  .   We  then  obtain 

2 
Lemma  2.   If  u  has  second  derivatives  in  L  then 

w  =  ru  satisfies 


2  r  r^u^  Ou  dv  =  2  /  w^Lw  dv  =  /  /  (w^+|  ^|  ^) t^ 


Here  (P  is  an  arbitrary  domain  in  t^r,^',^  space,  dv,  do- 
are  volume  and  surface  elements,  w  is  the  normal  derivative 
of  w  in  the  space  direction.  Also 
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a)   On  that  part  of  9  with  t  constant^  the  integrand 


is  w^  +  I Vw| ^. 


b)  On  that  part  of  (P  invariant  in  time  the  integrand 
is  -  2w^w^. 

c)  On  that  part  of  ^  that  is  characteristic  the 

1     ?         2     ? 
integrand  is  (w  +  (  |Vw|   -  w  ) ) ,  where  w  is  the 

^/2   ^  "^  ^ 

derivative  along  the  bicharacteristic . 

d)  Where  ff  is  space-like  the  integrand  is  positive 

definite  in  u  if  t  is  positive. 

n 

2.   The  Kelvin  Transformation. 


We  consider  the  Kelvin  transfonnation 


(1)   X'  =  x/(r^  -  t^),  t'  =  t/{r^  -  t^)  or  r'  =  r/r^  -  t^ 


leaving  c;',"*'  unchanged.   It  maps  the  exterior  of  the  cone 

r'^  =  t  into  itself  taking  the  origin  into  infinity  and 

2    2 
the  point  (^'^,0)  into  the  origin.   The  cone  r  =  t  is 

2   2      2    2  -1 
mapped  into  a  "cone  at  infinity"  since  r  -t  =  (r'  -t'  ) 

The  cones  r  +  t  =  k  are  mapped  into  the  cones  r'  t  t'  =  l/k. 

Thus  the  "cone  at  infinity"  orthogonal  to  r-t  =  k  is  mapped 

into  r'-t'  =  0.   To  find  how  the  operator  L  transforms  we 

note  that  A  is  invariant  and  b/ht  +   c)/or  =  (r'  r  t'  ) 
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{b/bV    :^   S/Sr').   Hence  h'^/ht^   -   S^/^r^  =  (r'^-t'^)^ 
(a^/St'^  -  b^/br^^).      Thus  rOu  =  Lw  =  (r'^  -  t'^)L'w  in 
obvious  notation.   Thus,  setting  u'  =  w/r'  we  see  that  if 
u  =  w/r  satisfies  Q  u  =  0  then  u'  satisifes*  O'u'  =  0. 

Futhermore,  Lenuaa  1  or  2,    the  energy  identity  can  be 

2    2 
applied  to  u'  in  the  primed  space  in  the  : mage  of  r  >  t  , 

♦       2     2 
i.e.  in  any  subdomain  of  r'   >  t'  ,  and  then  transformed 

back  into  the  unprimed  space  to  yield  an  identity  for 

r' 

u  =  —  u'  in  the  case  of  Lemma  1  or  directly'-  for  w  in  the 

case  of  Lemma  2.   However  it  is  awkward  to  transform  surface 
integrals  from  the  primed  to  the  unprimed  variables.  We 
note  instead  the  transformation  of  the  voliune  integral 

I  =  Tu^iQu'dv'  =  /  w^.Lw'  -^  =  /  (2rtw^,  ^-  (r^+t^)w^)Lw  r"^dv 
D'  <?'        ^'    <P 


2  A  2 tx  •  Vu  -:-  ( r ^+t  ^ )  u  +  2tu )  O  udv 


Furthermore 


2tx-Vu0u  -    (2rtu^(x-Vu))^  +  div   (tu^x)    +  2  div   (t(x-Vu)Vu) 


div   (t|Vu|^x) 


*See    L5J 
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is  a  quadratic  form  in  Vu,  u,  as  is 

(r^+t^)u^Ou  -  I  ((r^-i-t^)(u^+|Vu|2))^  +  div  ( (r2+t2)u^Vu) 
and  also 


2tul7u  -  2(tuu^)^  -  div  (2tuVu)  +  (u^)f 


In  general  a  quadratic  form  in  Vu^u,  cannot  be  a  divergence 
expression.   It  follows  therefore  that  when  the  three  above 
expressions  are  added  the  sum  must  vanish  since  we  knov;, 
from  the  primed  space,  that  the  sum  must  be  a  divergence. 

Thus  we  obtain 

2 
Lemma  3 .   If  u  has  second  derivatives  in  L  then 

I    (2tx-Vu  +  (r^+t^)u^+2tu) a  u  dv  =  T  //  2tu^x-Vu  +  i 

(r^+t^)u^+(Vu)2)  +  2uu^  -  u^  }  t^^  -  {   t(u^  -  (Vu)2)x.n 
+  [2t(x-Vu)  +  (lx|^+t^)u^H-2tu]u^}  x^  }  ds 


where  n  is  the  unit  space  normal  out  of  the  surface. 

It  is  also  clear  that  the  restriction  in  the  derivation 

2    2 
of  this  result  to  the  domain  r  >  t  may  be  dropped. 


.SJ)f. 
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We  can  obtain  an  identity  for  w  corresponding  to 
Lemma  2  by  adding  the  appropriate  divergence  expression  which 

in  this  case  Is  J    {-  |  (w^r'^lr^+t^)  )^^  +  -|  (w^r'^lr^+t^)  )^^> 


D 


dSdOdrdt  and  we  obtain 

2 

Lemma  4.   If  w  has  second  derivative   in  L  then 

J   (2rtw^+(r2+t2)w^)Lw  ^  =  J{l  ^   (r+t )2(w^+w^)2  +  ^  (r-t)^ 
(w^-w^)2  +1  (r2+t2){|vw|2  -  w^)  }  t^  -    ft(w^  -  |Vw|2) 
?.it  +  (2trw^  +  (r2-i-t2)w^)w^>  x^  }  -^| 


This  identity  was  derived  in  [2]. 

The  surface  integrand  is  positive  definite  for  space- 
like surfaces  with  t  >  0  and  involves  only  surface  deriva- 
tives on  characteristic  surfaces. 

The  last  statement  can  be  derived  by  noting  that  the 
integrands  of  the  surface  integrals  in  both  Lemmas  2  and  h 

involve  only  derivatives  of  v/.   Now  /  v;.  ,L'w  — ^  = 

(2rtw  +  (r  +t  )w. )Lw  ■^.   Hence  the  difference  between 

the  two  corresponding  surface  integrals  vanishes  identically, 
If  the  integral  over  S'  is  transformed  directly  into  an 
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integral  over  S  we  see  that  we  have  then  two  equal  integrals 
over  S  whose  integrands  are  quadratic  foiTTis  in  Vw  and  w,  . 
Hence  we  would  have  an  identically  vanishing  integral  over  S 
whose  integrand  is  of  the  form  Q-nx  +  Pt  where  Q,  P  are 
quadratic  forms  in  Vw^  w,  .   Thus  /  (div  Q  +  ^-t-)^^  vanishes 
identically.   But  since  Q  and  P  are  quadratic  forms  in  Vw, 
w.  it  follows  that  Q  =  0  and  P  =  0.   Hence  the  integrands  of 
the  two  surface  integrals  are  the  same.   Hence  the  surface 
integral  of  Lemma  4  is  what  vje  vjould  have  obtained  by  trans- 
forming directly  on  the  surface  integral  of  Lemma  2. 

If  the  surface  in  the  unprimed  space  is  characteristic 
the  image  in  the  primed  space  is  also  characteristic  and 
the  integrand  of  Lemma  2  involves  only  characteristic 
derivatives  in  the  primed  space  and  hence  after  transfor- 
mation only  characteristic  derivatives  in  the  unprimed 
space  as  we  v/anted  to  show. 

4.   Other  identities. 


It  might  now  appear  that  we  could  obtain  still  another 

identity  by  applying  the  Helmlioltz  transformation  T  to  the 

2 
primed  space.   However  T  is  the  identity  transformation 

and  hence  no  new  information  is  obtained. 

Every  invariant  transformation  leads  to  an  identity 

but  no  other  new  identity  can  be  derived  from  transforming 

the  classical  energy  identity.   For,  the  remaining  invariant 


■^-    srl-j   "ic    ■    ■■ 


-•i-OfiW 


rj 


V  C^  ■  \   . 
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transformations  are  a  translation,  rotation  or  stretching 
of  the  coordinates.   From  the  identity  in  Lemma  1  it  is 
clear  that  a  translation  produces  the  same  identity  and  from 
Lemma  2  a  rotation  does  the  same.   Stretching  the  variables 
x'  =  kx,    t'  =  kt  also  leaves  the  identity  of  Lemma  1  invariant, 

However  when  a  translation  or  rotaticn  is  applied  to 
the  coordinates  one  obtains  from  Lemma  K   nev;  identities. 
The  most  interesting  comes  by  taking  t'  =  t+c  in  the  identity 
of  Lemraa  k   and  equating  coefficients  of  c.   Thus  one  finds 

/jT  w(rw,,+tx^^)Lw  r"^dv  =  J  {  (^(r+t)  (w^+w^)^  +  llr-t)  (w^-w^)^ 

+  t(|Vw|2  -  w^J)t^  -  ((w^  -  IVwi^)  x.n 


+  2(rw^  +  tw, ))w^)xn  r   dS, 


This  identity  was  used  in  [l]. 

A  shift  of  the  origin  in  space  or  a  rotation  of  the 
axes  also  lead  to  ne^^r  identities.   Hov;ever  stretching  the 
variables  plainly  leaves  these  identities  invariant. 

5.   Hyperbolic  systems. 

A  somewhat  similar  principle  can  be  applied  to  a  vector 
solution  u  of  the  sj-mmetric  hyperbolic  system 


-lo- 


ft^ 
^t  =  ^  ^x 


X-,     X, 


X     1    2 
where  x  =  (x, ,  . .  .x  )  and  A  =  (A.   ^  A   ...)  are  constant 

matrices . 

Energy  conservation  yields,  over  any  domain  D, 

0=1      2u(u,  -  A^u  )dv  =  /  (u^t   -  u(A^-x  )u)do' 
J  t      X     J^n        n' 

D  D 


where  t  is  the  time  component  of  the  normal  and  x  is 
the  space  component. 

This  system  is  invariant  if  length  and  time  are  stretched. 
Hence  if  u(x,t)  is  a  solution,  so  is  ullacjkt)  and  hence 


Su 


=  X'Vu  +  tu,  is  also  a  solution. 
k=l  ^ 


Applying  the  identity  one  obtains: 


0  =  I    {   (x-Vu  +  tu^)^t^  -  ii?-^u  +  tu^)(A^x^)(x-Vu+tu^)  }   dcr 

From  this  one  may,  for  example,  conclude  that  if  u  has 
initially  compact  support  then  in  any  finite  region  the 

/  u,  |dx|  decays  like  — p  since  /   (x«Vu  +  tu,  )  d<^  is  bounded. 

'^  t=T 


:.  ^J 


f.i- 
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This  method  can  in  fact  be  used  on  the  wave  equation 
and  the  identity  of  Lemma  4  obtained  after  some  further 
integration  by  parts. 

II.  Application  of  the  Kelvin  transform  to  incoming 
and  outgoing  waves. 

In  [3>^]   an  outgoing  solution  u  of  the  v/ave  equation 
in  free  space  is  defined  as  one  X'/hich  vanishes  identically 
for  r  ^   t+k,  t  21  0*   By  the  Kelvin  transformation, 
r  =  r'/ 


=-  v-t/pi  _^\    ^      ^  -   t'/i"'  -t'   and  thus  the  corresponding 


j|x|l    =   /    ((j)     +    I V-^  I    )ldx|.      The   corresponding  Cauchy  data 


solution  u'  =  ^  to  an  outgoing  solution  vanishes  for 
r'  >  -t'+k"^,  t'  >  0. 

We  introduce  the  Hilbert  space  of  Cauchy  data  X  =  (v,  <{)) 
found  hy   the  closure  of  data  of  compact  support  in  tiie 
unprimed  space.   The  norm  is  the  energy  norm 

r 
J 

for  a  solution  w  of  Lx'j   =  0  will  be  (r.,  rtj))  and  the  identity 
of  Lemma  2  suggest  the  norm  ll}^!  ,  =  /  { (r({))^+|  V(rv, )  i  ^)  -^ 

which  we  call  the  v/-norm.   Since  the  space  of  data  is  found 
by  closing  the  space  of  data  of  compact  support  these  two 
norms  are  the  same . 

On  the  other  hand  the  Cauchy  data  ((j),v  )  a-re  carried  by 
the  Helmholtz  transformation  into  the  data 


A.   -      T 
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(-v',(l)')  with  V'  =  ^  V'   (-^,  -^),   <t)  =  -At  <!)  (-^,  -^) 

but  these  data  do  not  necessarily  form  an  element  of  the 
H' -space,  i.e.  Cauchy  data  of  finite  energy  nor^n  in  the 
primed  space. 

Hov;ever  there  do  exist  subclasses  of  data  in  the 
space  H' . 

It  is  also  clear  that  any  set  of  data  can  be  split 
into  its  outgoing  part  and  orthogonal  complement  by  making 
use  of  the  geometry  of  the  primed  space.   In  the  primed 
spc-ce  the  outgoing  wave  vanishes  on  the  cone  r'+t'  =  u 
and  the  orthogonal  complement  vanishes  in  r'  <  t' .   The 
two  components  of  the  initial  data  are  found  by  the  fol- 
lowing algoritlim.   First  solve  in  the  primed  space  the  Cauchy 
problem  for  the  corresponding  data  up  to  the  time  t'  =  k/2. 
Call  the  data  at  t'  =  k/2,  X^.   Then  solve  two  Cauchy 

problems  backv/ards  in  time  v;here  the  first  problem  has  the 

r' 
same  data  as  X,  for  IXJ  <  k/2  and  zero  outside  and  the  second 

problem  has  the  same  data  as  Xp  for  lX\'  ^   k/2  and  zero  inside, 

The  two  sets  of  corresponding  data  at  t  =  0  are  the  outgoing 

part  and  the  orthogonal  complement.   The  first  set  yield 

the  data  of  the  outgoing  component. 

To  make  the  argument  rigorous  we  note  the  follov;ing. 

lie  apply  the  identity  of  Lemma  4  in  the  primed  space  to 

the  slab  0  _;^  t'  <_  k/2  and  replace  w  by  aw^+bw^.   The 
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coefficient  of  ab  in  this  identity  gives  us  a  "scalar 
product"  identity  in  the  primed  space.   Now  if  w-.  and  Wp 
are  the  two  solutions  to  the  two  Cauchy  problems  described 
above  this  identity  reduces  because  of  the  choice  of  data 
on  t'  =  k/2  to 


r  2 

■-'  ('\.''w-j^-V'w2  +  w^^jWg^,  )dv' 

t"'=0 


If  this   is   ti^ansformed  to  the  unprimed   space  we  have 

r 
0  =   /      (Vw^-VWg+vr^^  ^gt^^"^ 

t=0 

or  the  two  sets  of  data  are  orthogonal  as  required. 

III.   Geometrical  optics  with  no  shadow. 

We  want  finally  as  an  application  to  sho^\r  that  geo- 
metrical optics  yields  asymptotic  solutions  to  the  reduced 
wave  equation.   The  problem  we  consider  is  the  following. 
Let  V  be  the  outgoing  solution  of 


AV  +  oo^V  =  5(x-a) 


which  vanishes  on  B  where  B  is  q)  star-shaped  with  respect 
to  a  point  inside  it  X'Jhich  we  take  as  origin.  We  "shall  show 
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that  V  is  given  by  the  Fourier  transform  /  ue'''"^  dt  where 
u  satisfies 


Du  =  0,  u  =  0  on  B,  u.  =  0,   u  =  ~  b{x-t) ,      for  t  =  0. 


Furthermore  this  transform  in  turn  is  asymptotic  to  the  formal 
expansion  in  en  found  on  taking  out  the  contributions  due 
to  the  singularities  of  u. 

Consider  the  solution  u  described  above.   It  is  given 
by  ico"  (p-t)&(p-'t')p~  ,  where  p  =  |x-S|,  in  the  region  bounded 
by  +  =  0,  the  body  cylinder  B  generated  by  the  body  in 
time  and  a  characteristic  surface ,  R,  formed  by  the  reflected 
rays  of  the  cone  p  =  t.  Across  this  characteristic  surface 
u  will  have  a  singularity  of  the  saine  type;  that  is,  if  D 
is  the  interior  of  the  reflection  surface  R,  then  in  the 
complement  of  D  exterior  to  B,  the  solution  is  given  by 


(1)      u  =  ioj"^  |{p-t)b(p-t)p"^+X^(il)5(U+X2S2(^)  +  ...  } 


*   * 
Here  X„,  X^  are  smooth  functions,  ^  is  the  normal  distance 

from  R  and  Sp  satisfies  /  ^'5(^')d^  =  o{>i)    etc.   For  details, 

see  [5]'   In  the  interior  of  D,  u  is  a  solution  of  the 

wave  equation  which  vanishes  on  B  and  is  a  smooth  function 


CD 


-1 


([)(^,t)  on  R.   Furthermore  on  R  the  integral 


-  15  - 


(2)   K{i»    =/ff  ^(r+t)2((r(t))^+(r(t))^)2  +  ^  (r-t)2(  (r(t))^-(r(t))^)' 
-    f  t((r(t))2  -  lV(r(l))l2)3f.r!  +  (2tr(r(t))^ 


is  bounded. 

Corresponding  statements  are  true  of  u,,u, ,,  etc. 

A  modification  of  Theorem  I  in  [2]  which  is  given  in 

r  2 
the  appendix  shovjs  that  cd  /  t  udt  converges  and  is  bounded 

T 

in  terms  of  K.   Here  T(x^y,  z)  is  the  value  of  t  vihere  the 
line  X  =  constant  cuts  the  reflection  surface  R.   Hence 


/       icDt . , 
/     ue   at 

T(x,y,z) 

CO 

exists.   Similarly  co  /  u5e~  dt  exists  as  do  the  integrals 


of  u, , .  And  all  these  integrals  are  bounded  independent  of 

CO. 

We  can  now  evaluate 


■^(,i<^'i) 


16  - 


00 


r    loot  T ,       .  r   -,  loot 

o)    /   ue        dt   =  -1    ;    ude 


T 


asymptotically  by  integrating  by  parts.   The  asymptotic 

icoT  /  *   —1  *    \ 
expansion  will  be  of  the  form  e    (X„+oo  X,  . .  . )  and  the 

remainder  will  be  of  order  co"  .   On  the  o-^her  hand 

T 
r     .   , 

o:>   I   ue   dt  has  a  similar  expansion  for  o  =  (x'-^|  =|=  0 
0 

CO 

ue'""'"  dt  is  asymptotically  equal  to  an 
0 

expression  of  the  form  e   '^(X^-Kd"  X,...)  where  X^^X  are 
functions  of  2. 

VJe  nov:  have 
Theorem:   If  u  is  a  weak  solution  of  0  u  =  0  exterior  to 
B  and  satisfies  u  =  0  on  B,  u  =  ioo"  6(x-a),  u,  =  0  for 

CO 

r   icut    exists  and 
t  =  0  then  /  ue   dt/is  asymptotically  of  the  form 
\J 
0 


e 


"^   '(X^-to"  X,...)  where  t(3?)  is  the  value  of  t  on  the 
characteristic  surface  formed  by  the  reflection  of  the  cone 
|x-1l|  =  t  on  the  body  B. 

We  must  finally  show  that  the  desired  function  V  is 


in  fact    /  ue""""^  dt.      VJe  note   that 


-  1?  - 


1^  oo 


0 


CO 

I 

0 


=  -  TAue-i^t--    "^" 


By  the  same  methods  as  in  [2]  and  noting  the  singularity  in 
u  at  t  =  t(x)  one  can  show 

00  00 

^Aue-i^^dt  =  A  j  ue-^^^dt. 
0  0 

Herce  the  integral  satisfies  the  reduced  wave  equation  with 
the  non-homogeneous  term  6{St-^) .      The  integral  vanishes  on 
the  body  B  since  u  =  0  there.   It  remains  to  show  that  the 
solution  is  outgoing,  i.e.  that  it  satisfies  the  Sommerfeld 
radiation  condition.   Again  consider 


CO 


/  ue    at  =  /  ue    dt  +  /  ue    ( 

0  0  T 


The  first  integral  consists  of  two  terms  e  '  /p  and  a 
term  of  the  form  e'^^'''X(T)  which  separately  satisfy  the 
Sommerfeld  radiation  condition  since  t  ~>  r  at  infinity  and 
X  becomes  independent  of  r.  The  second  integral  is  split 

*A  function  f  satisfies  the  Sommerfeld  radiation  condition 
if  r(f  -  icijf )  =  o(l)  as  r  — >  °«. 


.c"-} 


,  '.V      <~—    i    a .- 
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by  setting  u  =  ^j+u^  where  u^  is  the  solution  of  f?  u  =  0 
in  D  which  has  the  same  data  as  u  on  R  and  u^  is  the  solu- 
tion  given  by  the  retarded  potential  on  B,  i.e. 

Ug(^,t)  =  /  i  [  ^  ]dT  v;lth  R  the  distance  from  {x,t)   to 

the  variable  point  on  B  and  [   ]  means  the  retarded  value  on 
B   .   The  same  ai-gument,  slightly  modified,  as  that  in 
[2]  using  (vi)  of  the  modifieu  theorem  in  the  appendix 

shows  that  /  u^e"   dt  satisfies  the  radiation  condition. 
On  the  other  hand  /  u^e""^'""  dt  also  satisfies  this  condition 

T 

by  a  similar  argument  involving  the  representation  of  u^ 

by  mean  values  over  R. 

Hence  V  is  given  by  the  Fourier  transform  /  ue"'^''^  dt 

0 

and  thus  has  the  desired  as;^Tnptotic  expansion.   Substi- 
tution of  the  form  in  the  differential  equation  yields  the 
expansion  explicitly.    It  could  also  be  computed  directly. 


lO: 


■-.  1 
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Appendix . 

We  prove  here  a  modification  of  the  main  theorem  of 
[2]  that  is  necessary  for  the  results  of  Part  III. 
Theorem;   Let  -(>   be  a  strongly  star-shaped  infinite  body 

which  satisfies  /  — o  dor"  <  —  where  R  is  distance  from  x 

to  the  point  on  the  body-^  ,    S  is  the  minimiom  distance, 
£  is  a  constant,  k  is  a  constant,  and  dcr  is  surface  element 
on  Jf   .   Let  D  be  the  region  bounded  by  the  cylinder 
generated  by  Ms   in  time  and  a  characteristic  surface  R 
whose  generators  go  to  infinity  as  t  — >  ":   Let  u  be  a 
smooth  solution  of  0  u  =  0  in  D  which  satisfies  u  =  0  on 
and  u  =  (j)  on  R  where  (|)  satisfies 

]  {i\   (r+t)2((r(l))^+(r(l))^)2  +  |(r-t)2{  (r(l))^  -  (rct))^)^ 


|(r2+t2)(|V(r(l))i2  -  (r(t))^)  }   t^ 


-  f  t((r^)2  -  |V(r^)|2)^.^ 


+  (2tr(r^)^+(r2+t2)(r(t))^)(r^)^  1  ^n  ^Z"  ^ 


=   CO. 


K((t)) 


X  r- 


son  • 


3.1Gr,i^iC- 


f::.'i  iU'a 


--/       C*  •;  — -u  C  ■ 


-i      iJtw^, 


iU'      V- 


4     I  .;v  ^-^  ^ 
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Then  W  =  ru  satisfies  the  following  inequalities 


i)        /  -^  (  lVw|^  +  v/^)dv  <  ^  t^  for  t  large  enough, 
A  ^ 

where  A  is  a  fixed  region  in  space. 
Furthermore 


ii)  lu(x,t)i  -   ^ 

26^t 


<  ~— -  , 


r        r^    -,..  r.  K, 


iii)  J    t2  (|H)2^t  dc  <  ^  , 


where  a  depends  on  the  shape  of  the  body,  K,  on  K  and 
6  is  the  miniimun  distance  to  the  body^,  aiid 


iv)  rt\^(x,t)dt  <  ~  1+e 


00 


v)        r^""^  j  t2(u^.+u^)^dt  <  K^ 

0 

where  q  is  any  derivative  of  u. 

This  theorem  is  proved  by  proving  modified  forms  of 
Lemmas  1-6.   Lemma  1;,  [2],  is  another  formulation  of 
Lemma  k.     The  modification  is  to  consider  the  domain  D 
instead  of  a  slab  0  <^  t   <_  t^  .      In  Lemmas  2  and  3  the 
estimates  on  the  right  are  all  replaced  by  suitable  multiples 


•M,  f-  fJV 


e  ^: 


seji^i,-. 
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of  K.   In  Lemma  h,    Uj  is  a  free  space  solution  with  the 
same  data  on  R  instead  of  the  free  space  initial  value 

solution.   In  Lemma  5,  A/ 5  is  replaced  by  K/5  and  in 

-*    2  2+e 

Lemma  6  ,  r  becomes  r    and  A  becomes  k. 


^  2  4 

*In  [2 J  the  factor  r  should  read  r  .   In  the  proofs  of 

Lemmas  5  and  6  an  obvious  change  in  the  use  of  Schwarz' 

lemma  is  reouired. 


A. 
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